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ABSTRACT 


The  equations  of  motion  for  a  single  particle  traveling  in  a  constant 
density  atmosphere  are  derived.  The  aerodynamic  drag  on  the  particle 
and  the  atmosphere-particle  energy  transfer  resulting  in  loss  of  particle 
mass  are  considered  in  the  derivations.  It  is  assumed  as  an  initial 
condition  that  steady-state  ablation  is  occurring.  The  equations  are 
solved  using  a  Datatron  205  Digital  Computer,  and  the  solutions  are 
discussed  in  some  detail.  Emphasis  is  placed  on  determining  particle 
size  and  absolute  luminosity  from  measurements  of  distance  versus  time. 

Micron-size  particles,  which  travel  at  velocities  in  lower  meteor 
range  of  10  to  20  km/sec,  are  produced  by  impact  of  spherical  steel 
pellets  on  a  steel  target.  Preliminary  experiments  were  conducted 
using  these  particles  in  a  controlled  atmosphere.  The  particles,  being 
luminous,  were  detected  by  photomultiplier  tubes.  The  leading  edge  of 
a  cloud  of  particles  was  detected  and  velocities  up  to  15  km/sec  were 
measured.  By  applying  the  theory  to  deceleration  measurements,  the  size 
of  the  particles  was  estimated  at  approximately  1.0  micron  diameter. 

An  improved  vacuum  firing  range  was  designed  to  correct  for  the 
vacuum  and  size  limitations  of  the  original  system.  An  experiment  is 
proposed  to  utilize  the  improved  system  to  detect  and  measure  individual 
particles.  Data  from  the  experiments  can  be  compared  with  theory  and  the 
results  applied  directly  to  determine  in  detail  the  physical  phenomena 
occurring  in  meteor  flight. 
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I.  INTRODUCTION 


The  physical  theory  of  meteors  apparently  had  its  beginnings  in 

1922  and  1923  in  the  work  of  dpik  and  Lindemann  and  Dobson.  Since 

that  time  the  theory  has  been  expanded  and  considerable  observational 

data  have  been  collected.  Reports  of  original  work  and  summaries  of 
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the  work  of  others  are  given  by  Hoppe  ,  Whipple  ,  and  Hansen  .  The 

most  complete  and  the  most  recent  investigations  of  the  physics  of 
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meteor  flight  have  been  made  by  C)pik  and  Levin.  . 

The  question  naturally  arises  as  to  why  studies  are  being  made  of 
meteor  phenomena  and  why  they  are  important.  The  first  investigations 
were  carried  out  as  a  matter  of  curiosity  as  to  the  nature  of  meteors 
and  as  a  means  of  investigating  the  upper  atmosphere.  In  recent  years 
the  development  of  earth  satellites  and  space  vehicles  has  made  the 
presently  available  knowledge  of  meteors  extremely  important  and  is 
creating  an  ever  increasing  demand  for  further  insight  into  every  facet 
of  meteor  flight. 

Present  meteor  theory  is  based  on  assumptions  as  to  the  actual 
physical  processes  occurring.  Heat  transfer  to  a  vaporizing  surface, 
ablation  of  the  surface,  and  the  production  of  luminosity  are  not 
completely  understood.  As  a  result,  only  order-of-magnitude  estimates 
can  be  given  for  the  size  of  meteoroid  producing  a  given  visual  effect. 
Controlled  laboratory  experiments  are  a  necessity  before  the  problems 
of  meteor  physics  can  be  considered  solved. 


This  paper  presents  a  study  of  the  equations  of  motion  of  micron- 
size  particles  traveling  at  meteor  speeds  in  a  controlled  atmosphere. 

The  relations  between  size,  velocity,  distance  traversed,  and  luminosity 
of  the  particles  are  examined.  A  preliminary  experimental  investigation 

i 

is  reported  and  recommendations  are  made  for  future  studies.  Any 
information  gained  from  these  studies  can  be  applied  directly  to  actual 
meteor  problems. 

High  velocity  particles  used  in  the  experimental  work  are  produced 
by  shooting  a  projectile  from  a  gun  into  a  target.  As  a  result  of  the 
impact,  a  cloud  of  tiny  particles  is  produced  which  travels  away  from 
the  target  at  velocities  over  five  times  that  of  the  primary  projectile. 
These  particles,  henceforth  to  be  referred  to  as  spray  particles,  are 
very  similar  to  micro-meteors  in  both  size  and  velocity,  being  approx¬ 
imately  1/10  to  10  microns  in  diameter  and  traveling  between  5  and  20 
8  9 

km/sec.  ’  .  Spray  particles  can  thus  be  considered  to  be  experimental 

meteors  and  can  be  studied  as  such.  It  is  assumed  here  that  a  single 
particle  can  be  separated  from  the  cloud  and  studied  as  an  individual 
micrometeorite. 
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II.  THE  EQUATIONS  OF  MOTION  OF  AN  ABLATING  PARTICLE 


Derivation  of  Equations 

If  equations  are  to  be  derived  to  describe  a  problem  that  is  to  be 
tested  experimentally,  they  must  be  in  termsof  experimentally  measurable 
quantities.  In  spray  particle  studies,  time  of  flight  and  distance 
traveled  are  easily  measured;  therefore,  it  is  important  that  an  equation 
relating  these  two  parameters  be  determined. 

Let  it  be  assumed  as  a  first  consideration  that  the  high  velocity 

particle  is  not  losing  mass,  but  is  simply  being  decelerated  by  the 

atmospheric  resistance  to  its  motion.  With  the  particle  as  the  center 

of  the  frame  of  reference,  the  mass  of  air  swept  up  per  unit  time  is 

2 

pov,  resulting  in  a  time  rate  of  change  of  momentum  of  pov,  where  p  is 
the  density  of  the  atmosphere,  a  is  the  average  cross  sectional  area 
normal  to  the  direction  of  motion,  and  v  is  the  particle  velocity  measured 
in  the  laboratory  reference  frame.  The  decelerating  force  on  the  par¬ 
ticle,  m  can  be  related  to  the  change  in  momentum  by  the  factor  C^, 
the  drag  coefficient,  which  accounts  for  the  efficiency  of  the  momentum 
exchange  process.  The  drag  equation  is 


dv  ,  2 

m  dt  =  '1/2  CD  P°V 


Equation  (1)  can  be  integrated  to  obtain  velocity. 
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where  v  is  the  initial  velocity  of  the  particle  at  the  time  it  is 
o 

created.  The  distance  traveled  is  then 


vdt 
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Pv  t  +•  1 
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dt 


Perfifming  the  integration, 


x  = 


1 
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I  n  (1  +  (3v  t  ) . 
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(3) 


Assuming  that  the  particles  are  approximately  spherical  in  shape, 


pot’ 
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cmr  C 


3  =  - 
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2m 


2  (4/3t  8) 


3  °CD 
8  rS  ’ 


(4) 


where  8  is  the  density  of  the  particle  and  r  is  its  radius.  Equation 
(3)  generates  a  family  of  trajectory  curves,  so  that  by  changing  p  the 
various  members  of  the  family  are  obtained.  p  is  a  function  of  the 
four  parameters  p,  C  ,  8,  and  r.  For  experiments  in  a  controlled 
atmosphere,  p  and  6  are  known  leaving  only  C^./r  as  an  unknown.  For 
elastic  or  complet.ely-inelastic  collisions  between  meteor  and  air  molecules 
in  free  molecular  flow,  =  2 .  (A  more  complete  discussion  of  is 
given  later.)  The  curves  generated  by  Eq .  (3)  are  therefore  only  a 
function  of  particle  size.  An  experiment  can  be  performed  to  find  time 
of  flight  at  various  distances.  These  values  can  be  plotted  and  the 
resulting  curve  compared  with  Chose  of  £q ,  (3)  to  determine  particle  size. 
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Results  obtained  from  the  above  equations  would  be  In  error,  since 

the  particles  are  known  to  be  luminous  and  therefore  must  be  ablating. 

Consider  a  particle  traveling  with  a  velocity  v  through  a  medium  of 

constant  density  p.  Again  using  the  particle  as  the  center  of  the  reference 

frame,  the  mass  of  air  that  Is  intercepted  per  unit  time  is  opv,  where  a  is 

the  cross-sectional  area  normal  to  the  direction  of  motion.  The  kinetic 

tiE  3 

energy  swept  up  per  unit  time,  (^) ,  is  \  apv  .  If  a  fraction  y  of 
this  energy  goes  into  ablating  the  particle,  then 


and 


/i  3» 

=  y  a  apv  ) 


h 


dm  , 

dt 


dm  _  -  -g-  ycrpy' 

dt  h 


(5) 


where  m  is  the  instantaneous  mass  of  the  particle  and  h  is  the  heat  of 
ablation. 

The  drag  on  the  particle  is  given  by 


dv 

dt 


crpv 


The  deceleration  is 


dv 

dt 


j, 
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r  2 
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Time  can  be  eliminated  by  dividing  Eq.  (5)  by  Eq.  (6). 


or 


dm  =  ^mv 

dt  hCD 


dm  m  _ x. 
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h  C, 


vdv 
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If  the  particles  are  assumed  to  be  spherical,  which  is  a  good  supposition 
since  the  particles  are  small,  probably  molten,  and  more  than  likely 
rotating  which  tends  to  smooth  out  irregularities,  the  mass  of  a  particle 
of  density  6  is 

m  -  —Tfr35  , 


so 


dm 


m 


4»Tr2Sdr  3  dr 
4/3TTr35  r 


Substitution  in  Eq. (7) gives 


(8) 


dr  _  1  v  , 

r  TmT  vdv 


(9) 


Integration  of  both  sides  of  Eq.  (9)  results  in 


In  <~) 
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.  ,  2  2S 
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r  =  r  e 
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,  ,  2  2S 
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where  k  *  T  rrr-  flnd  r  and  v  are  initial  values.  Equation (6)  can  now 
o  00 

be  rewritten  in  terms  of  the  size  of  the  particle. 
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Substituting  the  value  of  r  from  Eq.  (10), 


dv 


Cnp  2  -k(v2  -  v2)  ,  0  -k(v2  -  v2) 
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where  a  =  — r  . 
o  r  o 
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It  will  be  later  found  that  the  luminosity  is  directly  proportioned  to  the 
time  rate  of  change  of  mass,  so  it  will  be  convenient  to  have  an  equation 
for  .  Again  assuming  a  spherical  particle 


4  3~ 

m  =  —  «r  o 


and 


dm  .  2C  dr 

-  =  4,tr  &  ~ 


(12.) 


From  Eq.  10 


dr  dv 

dE  =  2kvr  dF 


(13) 


Combining  Eqs.  (11),  (12),  and  (13) 


dm  Q  „  23 

dt  =  3,fk  CDP  r  V  • 


(14) 


With  the  previously  derived  equations,  the  luminosity  can  now  be  calculated 
in  terms  of  experimentally  measureable  quantities.  The  intensity  of 
luminous  radiation  is  usually  assumed  to  be  given  by 


_  t  dm  2 

1  =  ‘  2dt  v 


(15) 
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where  t  is  the  luminous  efficiency.  The  actual  processes  involved  in 
producing  the  visible  radiation  are  complicated  and  can  not  be  given 
accurately  by  the  simple  relationship  above;  however,  the  simple  equa¬ 
tion  seems  to  suffice  to  describe  the  essential  features  of  the  phenomenon,^ 
Observational  evidence  indicates  that  the  luminous  efficiency  will  not 
be  constant  but  will  vary  approximately  linearly  with  the  velocity'’’^’*'*' 
so  that 


Therefore , 


T  =  T  V 
O 


I 


o  dm 
2  dt 


3 

v 


(16) 


Substitution  of  Eq.  (14)  results  in 


I 


2*k  CD 


2  6 

pr  r  v 
o 


(17) 


Evaluation  of  Equation  Parameters 
Many  of  the  parameters  in  the  equations  of  motion  of  an  ablating 
particle  will  be  either  constant  or  restricted  to  a  certain  range  in 
any  experiments  that  can  be  performed.  It  is  thus  important  that  values 
or  ranges  of  values  be  determined  for  these  parameters. 

The  experiments  will  be  conducted  in  a  vacuum  tank  which  means 
that  the  atmospheric  density  inside  the  tank  can  be  controlled.  The 
range  of  values  for  p  will  thus  depend  on  the  pressures  that  can  be 
obtained  in  a  firing  range.  Most  meteor  phenomena  occur  in  regions  of 
the  atmosphere  where  the  mean  free  path  is  large  compared  to  the  meteor 
body.  For  the  easily  obtainable  pressure  range  of  0.1  mm  Hg.  to  10  mm  Hg. 
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the  mean  free  path  is  about  5  to  500  microns  which  is  of  the  same  magnitude  or 
larger  than  the  micron-sized  particles  that  will  be  studied.  Thus  for 


preliminary  investigations^ the  range  of  pressures  used  will  be  0.1  mm 

— 8  2/cm^ 

Hg.  to  10  mm  Hg.  or  in  terms  of  density  p  =  8  X  10  °  to 

1.5  X  10"5  g/cm  . 

To  simulate  iron-nickle  meteors,  the  first  experiments  may  be  conducted 

by  shooting  steel  balls  into  steel  targets.  The  density  of  the  spray 

3 

particles  will  therefore  be  that  of  steel,  so  8  =  7.8  g/cm  . 

8  9 

Previous  spray  particle  work  ’  has  shown  that  initial  velocities 
up  to  20  km/sec  are  present.  This  is  well  within  the  meteor  velocity 
range  of  11  to  73  km/sec.  The  range  of  values  that  will  be  used  in  the 
equations  is  10  to  20  km/sec.  Previous  work  has  also  shown  that  the 
particles  range  in  size  between  0.1  microns  £nd  10  microns  in  diameter. 

■S 

'  / 

These  values  will  be  used  for  r^  in  the  equations. 

The  coeficient  of  drag  to  be  used  in  the  equations  has  been  the 

12 

object  of  much  study.  Baker  concluded  that  the  molecules  that  are 

thermally  emitted  from  a  high  velocity  object  have  a  profound  effect 

,  ,  J  12 
upon  the  drag  and  cannot  legitimately  be  neglected.  Baker  and  Charwat 

found  that  by  incorporating  the  behavior  of  the  emitted  molecules  into 

the  analysis  of  a  sphere,  that  a  more  realistic  drag  coefficient  for  a 

sphere  has  the  form 

V 

CD  =  2  |  1  ~  DjB  +  (EQ  +  E^)  ,  (18) 

where  B  equals  the  ratio  of  the  diameter  of  the  sphere  to  the  mean  free 
path  of  the  emitted  molecules  with  respect  to  the  oncoming  molecules; 
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V  is  the  average  speed  of  the  emitted  molecules,  y  is  the  speed  of  the 
object;  =  0.23,  Eq  =  0.444,  and  =  0.65. 

Another  theory  claims  the  possibility  that  the  drag  coefficient 
is  a  function  of  velocity,  giving  for  the  drag  force. 


F0  " 


1  „  2- 
—  C_pv  a 
2  D 


[1  +  0.16s  t£“  ] 


(19) 


Here  v  is  the  velocity  of  the  high  speed  object,  y  is  the  heat  transfer 
coefficient,  and  s  the  object  velocity  in  units  of  10  km/sec.  Using 
this  equation  drag  is  increased  by  a  factor  of  one-half  for  a  velocity 
of  30  km/sec  and  is  doubled  at  60  km/sec. 

Cp,  the  drag  coefficient,  is  a  function  of  the  thickness  of  the 
shock  layer  or  cushion  of  air  that  appears  at  the  front  of  the  moving 
object,  and  also  of  the  shape  of  the  object.  All  sources  agree,  however, 
that  when  the  meteor  is  small  (in  the  micron  range)  the  air  cap  is  non¬ 
existent  or  negligible,  the  shape  is  unimportant  and  the  effects  of  the 
thermally  emitted  molecules  are  small.  Thus  the  drag  coefficient  takes 
on  a  value  of  approximately  two.  Since  the  particles  that  are  to  be 
studied  here  are  of  the  order  of  a  few  microns,  it  will  be  assumed  that 
Cp  is  exactly  equal  to  two. 

The  coefficient  of  heat  transfer,  y,  is  equal  to  the  coefficient 
of  accommodation,  K,  when  there  is  no  air  or  vapor  cap  to  shield  the 
high  velocity  particle  from  oncoming  air  molecules.  The  accommodation 
coefficient  is  equal  to  the  fraction  of  the  kinetic  energy  of  air  lost 
in  the  collision  process.  Let  it  be  assumed  that  v  is  the  velocity  of 
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the  oncoming  air  molecules  and  v‘  is  the  velocity  of  the.  reflected 
molecules  both  measured  relative  to  the  particle.  The  kinetic  energy 
lost  by  the  air  molecules  in  the  collision  is 


i  2  x  ,2 

2h.Y-_l-2H.Y_.. 

l  t2 
2PV’ 


1  -  (  — ) 
V 


(20) 


where  p  is  the  mass  of  a  molecule  of  air.  The  accommodation  coeficient 

is  therefore  only  dependent  on  the  incident  and  reflected  velocities  of 

the  air  molecules.  If  the  air  molecules  are  totally  reflected,  v'  =  -v 

and  K  *  0;  if  the  air  molecules  stick  to  the  particle,  v'  =  0  and 

£ 

K  =  1.  It  is  pointed  out  by  Opik  that  if  a  vapor  or  air  layer  is  present 

y  may  differ  considerably  from  K,  however,  for  this  study  of  micron- 

size  particles  the.  air  cap  is  negligible  and  y  ”  Since  it  is 

doubtful  that  the  air  molecules  will  be  completely  reflected  or  will 

stick,  a  raige  of  y  will  be  chosen  within  the  maximum  values,  y  will 

be  assumed  to  be  between  0.3  and  0.9.  A  theoretical  discussion  of  K 

£ 

can  be  found  in  “Opik  ,  and  a  discussion  based  on  experimental  investi¬ 
gation  of  K  is  given  by  Hansen^. 

The  heat  of  ablation,  h,  can  vary  over  a  wide  range  and  is  a  direct 

result  of  what  type  of  process  is  responsible  for  ablation.  For  a  high 

speed  particle,  ablation  may  be  by  vaporization,  fusion  and  spraying  of 

the  liquid,  sputtering,  fragmentation,  and  by  a  combination  of  any  or 

6 

all  of  these  processes. 

Sputtering,  which  unlike  vaporization  does  not  depend  on  temperature, 

is  a  process  whereby  the  impact  of  high-speed  atomic  particles  may  knock 

out  single  atoms  from  a  solid  or  liquid  surface.  Sputtering  seems  to  be 

14 

unimportant  except  for  very  high  velocity  particles  (50-80  km/sec) 
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The  heats  of  ablation  for  the  various  processes  can  be  estimated 

0*  £ 

by  using  the  following  values  for  iron  given  by  Opik  . 

Physical  Properties  of  Iron 


Melting  Point 
Boiling  Point 
Average  Specific  Heats 
Solid 
Liquid 
Gas 

Latent  Heat  of  Fusion 
Latent  Heat  of  Vaporization 
The  heat  of  ablation  for  melting  and 
temperature  of  283 J  K  is 


1800°  K 
3308°  K 

6 . 91  x  10^  erg/g°C 
£ 

6,66  x  10  erg/g°C 

Cv  =  2.23  x  106  erg/g°C 

Cp  =  3.72  x  10^  erg/g°C 
9 

2 . 69  x  10  erg/ g 
6  40  x.  10^  erg/g 

spraying,  h^,  assuming  an  initial 


h  =  6.91  x  106  (1800-283 )6C  +  2.69  x  109  =  (1.05  x  1010  + 

f  g  C  g 

0.269  x  10101  ^  =  1.30  x  1010 

g  g 

The  heat  of  ablation  for  vaporization,  h  ,  is 


h  =  h£  +  6.66  x.  106  (3508-1800) °G  +  6.40  x  1010  •^s=  8.84  x  1010  ■£££. 

v  f  g  C  g  g 


The  range  of  values  for  the  heat  of  ablation  is  therefore 


1.3  x  1010  ^  <  h  <  8.84  x  1010 

g  g 

From  the  above  discussion,  it  can  be  seen  that.  CL,  r,  h,  k,  and  a 

D  o 

will  be  approximately  constant  over  a  limited  range  of  experimental 
conditions.  In  integrating  the  equations  of  motion,  they  will  be  assumed 
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constant,  over  the  range  of  integration.  If  it  becomes  evident  that  the 
parameters  change  significantly  during  the  flight  of  one  particle,  a 
stepwise  integration  can  be  performed  to  approximate  the  actual  change 
in  conditions. 

Table  I  summarizes  the  parameter  values  that  have  been  chosen  as 
appropriate  for  the  present  study  of  micron-size  particles. 

Table  I  --  Values  of  parameters  to  be  used  in  equations  of  motion. 


Pressure 

(p) 

0.1  -  10  mm  H 

g 

Atm. Density 

(P) 

8xl0“8  -1.5xl0'5  g/cm3 

Particle  Density 

(S) 

7.8  g./cm3  < 

Initial  Radius 

<ro> 

0.05  -  5.0  micron 

Initial  Velocity 

<v 

10  -  20  km/sec 

Drag  Coefficient 

<v 

2.0 

Heat  Transfer  Coefficient 

(Y) 

0.3  -  0.9 

Heat  of  Ablation. 

CO 

1010  -  1011  erg/g 

k*  6K- 

.0025  -  .075  sec2 /km2 

V 

.  -2,  -1 

o  r  6 

o 

.41  x  10  -  77  m 

In  the  solutions  of  the 

equations 

of  motion  discussed  in 

section,  it  is  assumed  that  the  particles  have  been  heated  and  steady- 
state  ablation  conditions  attained  at  the  time  of  initial  observation 
of  the  flight.  A  more  complete  discussion  of  this  requirement  will  be 
given  in  a  later  report  along  with  extensive  calculations  of  trajectories 
under  a  wide  range  of  conditions. 
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Solution  of  Equations  and  Discussion  of  Results 


The  equations  of  motion  of  an  ablating  particle  may  be  integrated 

in  terms  of  exponential  integrals  to  obtain  particle  size  and  velocity 

as  a  function  of  distance.  This  is  discussed  in  Appendix  I. 

In  any  experimental  work,  it  is  convenient  to  measure  the  time 

taken  for  a  high  speed  particle  to  travel  a  given  distance,  hence 

Eq.  (11),  the  expression  for  the  particle’s  deceleration,  must  be 

integrated  to  arrive  at  an  expression  for  distance  as  a  function  of 

time.  Since  the  equation  is  not  directly  integrable,  and  numerical 

methods  must  be  employed,  it  was  decided  to  use  the  University  of  Utah 

Datatron  205  digital  computer  to  perform  the  integration.  The  values 

of  time,  velocity,  and  displacement  derived  from  the  deceleration  equation 

can  then  be  used  to  solve  the  remaining  equations  of  motion.  The  computer 

was  programmed  to  solve  for  the  following  quantities:  t,  v,  x,  r/r^, 

dm/dt  and  I/t  The  program  used  is  discussed  in  Appendix  II. 

The  manner  in  which  the  size  of  an  ablating  particle  changes  with 

time  is  shown  in  Fig.  1  which  is  a  plot  of  r/r  versus  time  for  various 

o 

values  of  the  size  parameter  a  *  C_p/r  6.  The  curves  are  drawn  for 

o  D  o 

v  *  15  km/sec,  y  =  0.6,  h  =  2  X  10^  joules/kg,  and  =  2.0,  thus 
2  2 

k  =  0.006  sec  /km  It  can  be  seen  that  as  the  value  of  aQ  becomes  larger, 

the  rapid  decrease  in  size  due  to  ablation  occurs  at  smaller  values  of  time 

and  therefore  also  at  smaller  values  of  displacement. 

It  is  interesting  to  note  that  allof  the  curves  of  Fig.  1  have 

identically  the  same  shape  and  all  reach  a  leveling-off  region  at  about 

r/r  =  0.25  where  ablation  continues  at  a  much  slower  rate.  If  a  value 
o 

of  r/r^  is  known  at  a  given  time  for  one  value  of  the  time  at 
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ratio  versus  time  for  an  ablating  particle  for  various  values  of  the  size 
leter  a  =  C  p/r  5.  v  =15  km/sec;  k  =  0.006  sec2/km^. 


which  r/ro  has  the  same  value  for  another  curve  is  found  by  multiplying 

the  given  time  by  the  ratio  of  the  second  value  of  a  to  the  first. 

This  allows  a  universal  curve  to  be  drawn  for  all  values  of  a  (v  and 

o  o 

k  are  constant  for  each  curve).  Three  of  these  curves  are  shown  in 

2  2 

Fig.  2  for  k  =  0.006  sec  /km  and  three  different  values  of  v  . 

O 

Fig.  3  shows  the  distance  traveled  by  an  ablating  particle  as  a 

2  2 

function  of  time.  Here  again  vq  =  15  km/sec  and  k  -  0.006  sec  /km  , 
and  the  curves  are  plotted  for  different  values  of  a  .  In  this  plot  also 
there  is  a  certain  regularity  in  the  shape  of  the  various  curves.  In 
this  case,  to  change  from  one  curve  to  another,  the  time  and  displace¬ 
ment  coordinates  are  multiplied  by  the  ratio  of  the  two  appropriate 

values  of  a  .  This  allows  a  universal  curve  to  be  drawn  for  all  values  of 
o 

ao  for  given  values  of  v^  and  k.  Three  such  curves  are  shown  in  Fig,  4  for 
different  values  of  vq.  For  a  given  value  of  a q  the  curve  deviates  very 
little  from  the  curve  of  constant  velocity  until  the  particle  radius 
has  been  reudced  to  approximately  70  per  cent  of  its  original  value.  At 
this  time  the  curve  bends  away  from  the  others  rapidly.  It  is  in  this 
range  that  the  size  of  a  particle  can  be  accurately  determined  experimentally. 

The  previous  discussion  has  been  limited  to  specific  values  of  k 
and  v  .  All  particles,  however,  will  not  have  identical  ablation  rates 
and  initial  velocities.  It  is  important  in  this  study  of  size  determin¬ 
ation  to  know  how  the  distance  versus  time  curves  are  changed  by  varying 
these  quantities.  Figure  5  shows  curves  for  specific  values  of  k  and 

a  for  three  initial  velocities.  It  is  clear  from  the  figure  that  for 
o 

small  values  of  time  each  curve  approaches  an  asymptote  dependent  on  the 


.  16  . 


al  curve  of  size  ratio  versus  time  for  an  ablating  particle  for  various 
of  initial  velocity,  k  =  0.006  sec^/km^. 


Distance  (cm) 


1  10  102  10 


Time  (psec) 

Fig.  3  -  Distance  versus  time  for  an  ablating  particle  for  various 

values  of  a  =  C  p/r  6„  v  “15  km/sec;  k  =  0.006  sec^/km2 
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(cm/m) 


19  - 


Distance  (cm) 


Fig.  5  -  Variance  in  distance  versus  time  curves  resulting  from 
a  change  in  initial  veljcitv. 


*  1.0  m  ,  k  ~  0. .'06  Fec^/km^. 


value  of  v  .  The  asymptotes  are,  of  course,  just  the  curves  for  constant 
o 

velocities  equal  to  v  .  Experimentally  it  should  not  be  difficult  to 
determine  the  initial  velocity  for  a  given  particle  since  the  deviation 
from  the  asymptote  in  each  case  occurs  in  approximately  the  same  distance 
range  but  at  widely  different  times.  An  ideal  method  of  determining 
initial  velocity  would  be  to  have  a  detector  placed  very  close  to  the 
target  to  measure  the  velocity  at  this  point. 

Figures  6,  7,  and  8  show  the  variance  in  a  distance  versus 
time  curve  for  various  values  of  k  for  initial  velocities  v  ,  of!  10, 


°  / 

15  and  20  km/sec.  The  difference  in  k  arises  from  different  hjfeat 

transfer  coefficients  and  heats  of  ablation.  The  non-ablation  case  of 

Eq.  3  is  represented  by  k  *  0.  Ablation  by  vaporization  is  shown  by 

the  curves  for  k  ■  0.0025  and  k  =  0.006,  and  ablation  by  melting  and 

spraying  by  k  *  0.05  in  Fig.  6.  Since  the  process  of  ablation  by 

spraying  predominantly  occurs  with  larger  particles,  values  of  k  over 

2 

approximately  0.01  (sec /km)  will  not  be  important  with  the  tiny 
particles  that  are  of  interest  in  this  study. 

In  the  range  of  k  used  here,  there  is  a  very  little  difference  between 
the  ablating  particle  and  the  non-ablating  particle  for  low  initial 
velocities.  Depending  on  the  accuracy  of  experimental  work,  it  may  be 
perfectly  all  right  to  determine  initial  size  from  the  equation  for  a 
non-ablating  particle.  However,  when  size  at  a  given  time  is  required, 
the  ablation  case  must  be  used.  For  high  velocities  k  becomes  more 
important  as  shown  by  Fig.  7  and  8.  In  the  following  discussion  it  will 
be  seen  that  the  value  of  k  is  very  important  in  the  study  of  luminosity 
and  will  have  to  be  considered. 
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Distance  (an) 


i 


) 


Time  (fisec) 


Fig.  6  -  Distance  versus  time  curves  for  an  ablating  particle  for 
various  values  of  k.  a  -  1.0  m"^;  v  *  10  km/sec. 


22 


Distance  (cm) 


Distance  (cm) 


Fig.  8  *  Distance  versus  time  for  an  ablating  particle  for  various 
values  of  k.  olq  =  1.0  m-1;  =  20  km/sec. 
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The  luminosity  of  a  high-speed  particle  is  given  by  Eq.  (17)  as 


I  =  3  / 2rtkCjj  pTQr2  (17) 

It  was  found  in  programming  the  computer  that  the  luminous  intensity 
could  be  obtained  more  conveniently  from  the  expression 


t  C  pr 
o  0  o 


=  3/2*k  (r/r  )2  V& 

O' 


(20) 


Equation  20,  differing  only  by  a  constant  from  Eq.  (17),  is  plotted  in 

Fig.  9  for  k  =  0.006  (sec/km)  and  v^  =  15  km/sec.  The  curves  show  a  marked 

difference  for  various  values  of  the  size  parameter  a  ;  thus,  the  luminous 

o 

intensity  of  a  large  particle  remains  constant  at  its  initial  value  for  a 

much  longer  time  than  for  a  small  particle. 

As  with  the  curves  of  Figs.  1  and  3,  a  universal  intensity  curve 

can  be  drawn.  The  same  intensity  for  one  curve  is  obtained  on  another 

curve  by  dividing  the  value  of  a ^  by  a  given  amount  and  multiplying  the 

time  by  the  same  amount.  Thus  one  curve  can  be.  drawn  for  specified  values 

of  k.  and  v  .  The  intensity  curve  for  any  size  particle  can  then  be 

found  by  multiplying  the  time  scale  by  the  required  correction  factor. 

Universal  curves  are  shown  in  Fig.  10  for  different  velocities. 

Figures  11  and  12  show  how  the  luminous  intensity  varies  with  changes 

in  v  and  k.  For  a  given  size  particle  there  is  a  difference  in  intensity 
o 

of  approximately  a  factor  of  10  when  the  initial  velocity  is  different  by 
a  factor  of  5.  There  is  also  a.  progressive  increase  in  the  time  rate 
of  change  of  intensity  as  initial  velocity  is  increased.  It  can  be  seen 
from  Fig.  11  that  for  intensity  studies,  the  initial  velocity  must  be  measured 
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Fig,  10  -  Universal  luminous  intensity  curves  for  an  ablating 
particle  for  various  valuer  of  initial  velocity, 
k  1=1  0,006  sec^/fo?/-. 


(km  / sec  ) 


0.1  1.0  10  102 


Time  (|j,sec) 

Fig.  11  -  Luminous  intensity  versus,  time  for  an  ablating  particle 
for  various  values  of  initial  velocity,  a  -  1.0  m"^-; 
k  =  0.006  see2/'km2.  ° 
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(km  /sec 


Fie  12  -  Luminous  intensity  versus  time  for  an  ablating  particle  for 

various  values  of  k.  Oq  »  1.0  aT1;  VQ  -  10  km/sec 
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experimentally  to  within  at  least  1  km/sec  before  any  accurate  conclusions  can 
be  drawn  on  luminosity. 

The  intensity  varies  with  k.  The  value  of  k  is  not  easily  measured 

by  experiment  and  so  must  be  calculated  on  the  basis  of  present  knowledge. 

Assuming  the  particles  are  ablating  by  vaporization,  k.  will  have  a  value 
9  2 

close  to  0.006  sec“/km  .  The  error  in  luminosity  using  k  =  0,006  when 

the  actual  value  is  0,002,5  is  a  factor  of  about  2.3  based  on  the  curves  of 

Fig.  12.  If  this  amount  of  inaccuracy  can  be  tolerated  a  value  of 
2  2 

k  •»  0.006  sec  /km  can  always  be  used. 

It  is  possible  that  the  size  of  high  velocity  particle  size,  can  be 
measured  by  an  experiment  that -measures  their  luminosity.  Since  the  luminous 
intensity  varies  strongly  with  particle  size,  an  accurate  measure  of 
intensity  versus  time  could  easily  indicate  size. 
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III.  EXPERIMENTAL  STUDIES  OF  SPRAY  PARTICLES 


High-velocity  spray  particles  produced  by  the  impact  of  a  projectile 
on  a  target,  provide  a  source  of  meteor- like  particles  and  allow  the  labor¬ 
atory  investigation  of  problems  of  meteor  physics.  With  these  particles, 
the  equations  derived  in  the  previous  section  may  be  tested  and  the  para¬ 
meters  evaluated. 

Earlier  work  on  spray  particles  by  J.  S.  Clark,  R.  R.  Kadesch, 

R.  W.  Grow,  W.  H.  Clark,  R.  E.  Blake,  and  E.  P.  Palmer  8>9>]-5>16  has  re_ 
suited  in  the  determination  of  many  of  the  properties  of  the  particles  and 
the  development  of  experimental  techniques.  In  particular,  it  has  been 
determined  that  the  particles  range  in  size  from  1/10  to  10  microns  and 
velocities  up  to  20  km/sec  exist.  Unpublished  work  by  R.  R.  Kadesch  in¬ 
dicates  that  higher  velocities  may  be  present  and  that  the  smaller  particles 
have  the  higher  velocities.  Limitations  in  the  experimental  equipment 
used  in  the  earlier  work  has  limited  the  investigation  to  the  size  and 
velocity  ranges  indicated.  In  all  earlier  spray  particle  work,  the  cloud 
of  particles  produced  in  an  impact  has  been  observed  as  a  whole.  Velocity 
measurements  have  been  made  of  the  leading  edge  of  a  moving  cloud. 

The  experimental  work  reported  here  has  been  aimed  at  obtaining  in¬ 
dividual  spray  particles  with  known  size  and  velocity  in  order  to  more 
effectively  investigate  meteor  phenomena  aid  hyperveiocity- impact  phenomena. 
Preliminary  experimental  work  will  be.  discussed  and  the  limitations  imposed 
by  the  equipment  used  will  be  considered.  An  improved  firing  range,  which 
has  been  built  will  be  described  and  proposed  experiments  in  the  new  range 

\ 
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will  be  presented. 


Design  of  Experiment  to  Measure  Particle  Size 

In  Section  II,  the  equations  of  motion  of  an  ablating  particle  were 
derived  and  solved.  The  next  step  is  to  obtain  experimental  data  to  com¬ 
pare  with  the  theoretical  work. 

Spray  particles  are  luminous;  therefore,  a  light  sensing  device 
can  be  used  to  detect  the  particles  as  they  pass,  and  a  series  of  these 
detectors,  together  with  appropriate  recording  equipment,  can  give  time 
and  distance  information  for  various  points  along  the  line  of  flight. 

Using  the  experimental  values  of  distance  and  time  together  with  the 
theoretical  work  of  the  previous  section,  particle  size  and  luminosity 
can  be  determined.  The  accuracy  of  the  measurement  depends  upon  the  time 
and  space  resolution  of  the  particle  detectors. 

In  experiments  with  spray  particles,  the  density  of  the  ambient 
atmosphere  becomes  very  important.  Very  small  particles  are  decelerated 
extremely  fast  in  a  dense  atmosphere  and  may  travel  less  than  a  centimeter, 
whereas  large  particles  in  a  rarefied  atmosphere  may  not  lose  enough  speed 
in  a  few  meters  to  make  deceleration  measurements  practical.  Thus  it  is 
essential  that  the  experimental  work  be  carried  out  in  a  controllable 
atmosphere. 

In  a  practical- sized  laboratory  vacuum  tank,  measurements  can  be 
made  over,  a  range  of  about  4  to  200  cm.  Fig.  3  shows  that  over  this 
range,  size  determinations  can  be  made  if  =  C^p/r^b  varies  between 
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about  1.0  and  10  m  .  To  study  a  given  size  of  particle,  the  value  of 
density,  p, must  then  be  controlled  such  that  0Lq  falls  between  these 
values . 

Using  this  method  of  determining  size  and  having  a  liir.ited  number 
of  detectors,  it  is  only  possible  to  determine  if  a  particle  is  or  is  not 
a  certain  size.  As  an  example  for  the  v^  and  k  values  of  Fig.  3,  a  par¬ 
ticle  whose  size  makes  the  quantity  0Lq  equal  to  10  m  ^  can  be  measured  if 
several  detectors  are  placed  in  the  distance  range  from  about  5  cm  to 

15  cm  from  the  target  area.  If  a  =  1.0  m  \  the  detection  equipment 

o 

would  have  to  be  placed  about  50  to  150  cm  from  the  target;  if 

a  =  100  m  \  the  range  would  be  0.05  to  0.15  cm.  Thus,  for  a  given 
o 

positioning  of  the  detection  apparatus,  the  following  information  on 
particle  size  will  be  available:  The  particle  is  in  a  certain  size  range 
or  the  particle  is  larger  than  this  size.  It  cannot  be  smaller,  since  a 
smaller  particle  will  be  decelerated  rapidly  and  will  not  reach  the  detection 
equipment  with  sufficient  velocity  to  be  luminous. 

Experimental  Equipment 


Vacuum  Firing  Range 

Preliminary  experimental  studies  of  spray-particle  dynamics  were 
carried  out  in  the  vacuum  firing  range  shown  in  Fig.  13.  The  tank,  is 
6x1x1  ft.  and  is  made  of  1/2  in.  steel  plate.  There  are  two  large 
hatches  on  the  top  of  the  tank  (shown  by  dotted  lines)  and  five,  small 
openings,  one  at  each  end  and  three  along  the  front  face.  Two  3/4  in. 
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Fig.  13  Diagram  of  Vacuum  Tank  and  Firing  Range 


aluminum  plates  are  used  for  coverings  for  the  hatches  and  are  held  down 
by  screw  clamps.  Four  of  the  small  openings  are  covered  with  plexiglass 

for  viewing  purposes,  and  the  fifth  or  center  opening  is  connected 

through  a  blast  tank  to  the  gun. 

The  vacuum  pump  is  a  Welch.  Model  1403 -B  high  vacuum  pump,  rated  at 
0,00.5  mm  Hg  ultimate,  pressure  and  a  pumping  speed  of  100  liters  per 

minute.  It  is  connected  to  one  end  of  the  tank  by  a  3/4  in.  metal  pipe. 

The  pressure,  measuring  apparatus  consists  of  an  oil  manometer 
which  reads  down  to  approximately  10  cm  Hg  and  a  small  mercury  manometer 
reading  down  to  0.5  mm. 

The  gun  used  to  accelerate  the  primary  projectiles  is  mounted  on 
a  sturdy  wooden  table  15  feet  from  the  main  tank.  The  gun  barrel  consists 
of  two  24  in.  sections  of  1-1/4  in.  cylindrical  steel  stock  containing  a 
22  caliber  bore.  One  end  of  the  barrel  is  chambered,  for  a  220  swift 
cartridge.  The  gun  is  fired  by  means  of  a  solenoid  which  forces  the 
firing  pin  into  the  primer  at  the  rear  of  the  cartridge, 

A  blast,  tank  is  located  at  the  forward  end  of  the  gun.  It  con¬ 
sists  of  a  series  of  steel  plates  each  containing  a  hole  just  slightly 
larger  than  the  size  of  the  pellet  that  is  fired.  These  plates  act  as 
baffles  which  allow  the  pellet  to  pass  unhindered,  but.  show  a  high  re¬ 
sistance  to  the  powder  gases  following  the  pellet.  The  gases  thus  take 
longer  to  reach  the  ta.rget  area  tba.n  the  pellet  and  do  not  interfere  with 
the  experiment, 

A  4-1/2  in.  diameter  hollow  pipe  connects  the  blast  tank  to  the 
center  small  opening  on  the  front  face  of  the  main  tank.  This  pipe 
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makes  it  possible  to  obtain  a  vacuum  throughout  the  system,  including 
the  barrel  of  the  gun,  and  also  provides  a  shield  in  case  the  pellet 
does  not  travel  straight  down  the  range  but  careens  off  to  one  side. 

Electrical  connections  can  be.  made  to  equipment  inside  the  tank 
by  means  of  several  coaxial  connectors  located  at  various  positions  on 
the  top  and  a  few  general  purpose  connectors  on  the  front  face. 

Spray  Particle  Detection  Equipment 

The  luminosity  of  spray  particles  in  an  atmosphere  makes  it 
possible  to  use  light-sensing  devices  for  particle  detection.  The 
present  equipment  uses  three  photomultiplier  tubes  and  associated  cir¬ 
cuitry  to  detect  a  high-speed  particle  as  it  passes  and  to  give  an  out¬ 
put  proportional  to  the  luminous  intensity.  Figure  14  shows  a  circuit 
diagram  for  one  of  the  three  identical  ph.otomultiplier-tube  detectors. 

The  photomultiplier  tube  is  an  RCA  Type  931-A  which  has  a  maximum  current 
amplification  of  800,000  using  a  voltage  difference  of  1000  volts  d-c 
between  anode  and  cathode  and  100  volts  between  dynodes.  An  additional 
gain  of  3  is  proveded  by  a  one-stage  amplifier.  A  cathode  follower 
matches  the  output  impedance  to  that  of  the  coaxial  cable  used  to  transmit 
the  detection  information  to  equipment  outside  the  vacuum  tank. 

Three  photomultiplier  tubes  are  mounted  in  a  straight  line  and 
placed  in  a  rectangular  aluminum  box.  The  box  serves  as  a  light  shield 
and  holdsthe  detectors  securely  positioned. 

The  distances  between  photomultiplier  tubes,  and  from  the  first 
one  to  the  front  end  of  the  box,  are  known  to  within  0,1  mm.  At  the 
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Fig.  14  Circuit  diagram  of  spray  particle  detector 


front  end  of  the  box  a  double  wall  serves  as  a  shield  to  limit  the 
number  of  particles  entering  the  box.  The  wall  sections  are  3/8  in. 
apart  and  each  contains  a  small  hole.  Spray  particles  traveling  in  a 
horizontal  plane  and  in  a  position  to  pass  the  photomultiplier  tubes 
will  pass  through  the  double  hole  and  be  detected. 

It  is  desired  that  a  given  detector  responds  only  when  a  particle 
passes  directly  in  front  of  it.  To  restrict  the  field  of  view,  a 
cylindrical  aluminum  cover  shield  containing  a  1/4  inch  slot  is  placed 
over  the  photomultiplier.  A  plane  shield  with  an  adjustable  slit  is 
placed  immediately  in  front  of  the  slot  in  the  cover.  With  the  variable 
slit  adjusted  to  0. 5  mm  and  positioned  1.0  cm  from  the  face  of  the  photo 
multiplier,  the  viewing  angle  is  approximately  3°.  With  this  small  angle 
of  view,  the  particle  will  be  detected  as  it  travels  a  distance  of  0.6  mm. 
The  ideal  angle  of  view  is,  of  course,  0°,  but  as  the  slit  is  made 
narrower,  light  intensity  is  reduced  and  detection  becomes  difficult. 
Figure  15  is  a  photograph  showing  the  particle-detection  equipment. 

A  block  diagram  of  the  experimental  arrangement  can  be  seen  in 
Fig.  16.  Each  of  the  three  photomultiplier  circuits  is  connected  by 
means  of  coaxial  cable  to  the  vertical  input  of  an  oscilloscope.  A 
phototube  mounted  above  the  target  area  views  the  flash  occurring  at 
impact  and  sends  a  triggering  pulse  simultaneously  to  the  three  oscillo¬ 
scopes.  Outputs  of  the  detectors  are  displayed  on  the  oscilloscope 
screens  and  photographed  with  a  Dumont  oscilloscope  camera.  Three  points 
on  a  distance-versus-time  plot  for  the  particles  are  then  known. 
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(b) 

Fig.  15  Photographs  showing  spray  particle  detection  equipment. 

Photomultiplier  tubes  and  mounting  board  (a)  outside  of 
box;  (b)  in  position  inside  box. 
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Fig.  16  Block  diagram  of  spray  particle  detection  equipment 


Preliminary  Experiments  and  Results 

Preliminary  experimental  work  was  done  using  the  vacuum  system 

and  firing  range  described.  The  detector  box  was  positioned  in  the 

vacuum  tank  as  shown  in  Fig.  13.  The  tests  were  made  with  7/32  in. 

steel  balls  shot  into  a  flat  steel  target.  The  detector  box  was  placed 

so  that  the  input  holes  were  lined  up  approximately  in  the  plane  of  the 

target  and  at  right  angles  to  the  path  of  the  primary  projectile.  It 

9 

was  found  in  earlier  particle  work  that  the  highest  velocity  spray  comes 
off  at  very  shallow  angles  to  the  plane  of  the  target. 

The  first  successful  operation  of  the  complete  detection  system 
resulted  in  the  data  shown  below.  The  pressure  in  the  tank  was  1.5  cm  Hg. 


Shot  1 

Shot  2 

Shot  3 

X 

t 

t 

t 

Photomultiplier 

cm 

(u  sec) 

(u  sec) 

(u  sec) 

1 

12.8 

15 

17 

13 

2 

18.5 

25 

36 

25 

3 

27.1 

37 

>50 

38 

Shot  1 

Shot  2 

Shot  3 

Average  velocity 

between  origin  and 

tube  1 

8 . 5  km/ s  ec 

7.5  km/ sec 

9.8  km/sec 

Average  velocity 

between  tubes  1  and 

l  2 

5.7 

2.8 

4.3 

Average  velocity 

between  tubes  2  and 

l  3 

4.7 

<4.0 

4.3 

Oscilloscope  traces  of  a  typical  shot  are  shown  below.  The  data 
above  are  actually  for  a  cloud  of  spray  particles  and  not  for  a  single 
particle.  This  is  so  because  the  input  holes  were  1/4  in.  in  diameter 
and  permitted  many  particles  to  enter. 
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Traces  are  from  right  to  left. 


Output  1  Output  2  Output  3 

2  volts/cm  5  fisec./cm 

For  a  possible  estimate  of  the  size  and  velocity  of  the  particles, 
the  results  of  shot  3  were  compared  with  the  theory  presented  in 
Section  II.  From  Fig.  4,  it  was  found  that  the  three  experimentally  de¬ 
termined  points  lie  approximately  on  the  curve  for  =  15  km/sec  and 

-1  3 

aQ  =  C^p/r^  *  10  m  .  With  equal  to  2.0,  6  equal  to  7.8  g/cm  ,  and 

-5  3 

p  equal  to  2.4  x  10  g/cm  which  corresponds  to  a  pressure  of  1.5  cm  Hg, 

*  -  6 

rQ  is  calculated  to  be  0.62  x  10  m.  This  result,  although  not  exact 

since  the  data  are  for  a  cloud  of  particles  and  not  an  individual  one, 

8  9 

agrees  with  the  size  and  velocity  estimates  of  earlier  work.  ’  In  a 
cloud  of  particles  the  large,  initially  slower  particles  overtake  and 
pass  the  small,  initially  faster  ones.  Over  long  distances  this  would  have 
to  be  considered;  however,  over  the  range  of  detection  for  the  data  listed 
above,  the  effect  is  probably  unimportant.  Thus,  the  leading  edge  of  the 
cloud  is  made  up  of  particles  of  nearly  uniform  size  and  velocity.  From 
Fig.  2,  it  is  found  that  particles  with  an  initial  velocity  of  15  km/sec 
are  reduced  in  size  to  32%,  28%,  and  26%  of  their  original  radius  at 
the  time  they  pass  the  first,  second,  and  third  photomultiplier  tubes, 
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respectively.  Little  can  be  determined  concerning  luminosity  because 
of  the  unknown  number  of  particles  that  were  detected. 

The  input  holes  were  then  made  smaller  to  obtain  fewer  particles. 

An  estimate  of  the  actual  number  of  particles  being  detected  was  made  by 
using  a  microscope  to  count  the  number  of  craters  produced  in  a  highly 
polished  target  placed  directly  following  the  third  photomultiplier 
tube.  With  2  mm  diameter  shield  holes,  about  .50  craters  were  observed. 
The  fact  that  a  group  of  50  particles  can  be  easily  seen  is  of  great 
interest  in  indicating  that  a  single  particle  could  be  detected  under 
similar  conditions  with  only  slight  improvements  in  the  light  gathering 
capabilities  of  the  detectors. 

At  this  stage  of  the  experimental  work  the  old  vacuum  tank  was 
abandoned  and  work  begun  on  a  new  system. 

Design  of  a  New  Vacuum  Firing  Range 

With  experience,  it  became  evident  that  the  vacuum  range  used  in 
the  first  experiments  was  inadequate  for  single-particle  studies.  The 
main  difficulty  was  that  the  pressure  in  the  system  could  not  be  re¬ 
duced  enough  to  allow  small,  fast  particles  to  reach  the  detectors.  A 
second  difficulty  was  that  limited  space  prevented  the  most  advantageous 
placement  of  detectors  and  alignment  wa#  difficult.  To  overcome  these 
limitations,  a  new  vacuum  range  was  designed  and  is  now  nearing  com¬ 
pletion. 
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A  top  view  of  the  new  system  appears  in  Fig.  17.  The  main  tank 
is  4  ft.  in  diameter  and  23  in.  deep.  The  lid  is  a  large  circular  piece 
of  steel  plate  braced  with  crossed  I-beams  on  the  top.  It  is  positioned 
on  the  tank  by  means  of  a  counterweight  and  pulley  system.  Seven  ports 
around  the  sides  of  the  tank  can  be  used  for  various  purposes.  A  series 
of  5  ft.  lengths  of  10  in.  diameter  iron  pipes  is  connected  to  one  of 
the  ports.  The  gun  described  in  the  discussion  of  the  old  system  fires 
through  these  large-diameter  pipes  and  into  thelarge  main  tank.  An 
aluminum  plate  containing  a  1  in.  hole  at  its  center  is  placed  between 
each  section  of  pipe.  This  series  of  plates  acts  as  a  baffle  to  slow 
the  gases  that  follow  the  pellet  down  the  range.  Electrical  connections 
can  be  made  to  a  number  of  hermetically  sealed  connectors  mounted  on 
an  aluminum  plate  bolted  to  one  of  the  large  ports. 

The  vacuum  equipment  consists  of  a.  Stokes  Model  212H-10  vacuum 
pump  powered  by  a  5  hp,  220  volt-3  phase  motor.  The  pump  is  rated  at 
0.010  mm  Hg  ultimate  pressure,  with  a  pumping  speed  of  140  cubic  feet 
per  minute.  The  vacuum  connection  to  the  tank  is  made  with  3  inch  copper 
tubing  through  a  3  inch  gate  valve.  Pressure  in  the  tank  is  measured 
with  a  Kinney  Type  T.D.-l  McCloud  vacuum  gauge. 

The  tank  has  adequate  space  for  positioning  detection  equipment 
to  measure  particles  leaving  the  target  at  any  angle.  Extension  tubes 
can  be  placed  on  the  large  ports  for  long-range  measurements.  Tests 
show  that  pressures  down  to  0.1  mm  Hg  are  easily  obtained  which  should 
allow  the  observation  of  small  particles  traveling  at  higher  velocities 
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Fig.  17  Diagram  of  new  vacuum  tank  and  firing  range 


than  in  the  first  system. 


A  Proposed  Experiment  in  the  New  Vacuum  Range 
With  the  theoretical  background  of  Section  II  and  the  experience 
gained  in  the  first  experiments,  a  proposal  for  an  experiment  to  utilize 
the  capabilities  of  the  new  vacuum  firing  range  can  be  made. 

The  first  efforts  will  have  to  be  directed  toward  obtaining  a 
single  spray  particle,  since  the  theoretical  studies  were  made  for  one 
particle.  It  appears  that  the  simplest  way  to  do  this  is  to  reduce  the 
size  of  the  input  holes  to  the  detector.  As  the  holes  are  made  smaller, 
the  problem  of  alignment  of  the  detection  equipment  will  be  increased 
and  the  probability  of  getting  a  suitable  particle  will  be  decreased. 
Construction  of  a  precision  positioning  device  may  help  to  solve  this 
problem.  It  will  be  assumed  for  the  discussion  to  follow  that  it  is 
possible  to  obtain  a  single  particle. 

If  the  old  detector  box  with  its  fixed  spacing  is  used,  it  should 
be  placed  so  that  the  first  photomultiplier  tube  is  approximately  5  cm 
from  the  target.  Referring  back  to  Fig.  2,  it  can  be  seen  that  5  cm 
to  15  cm  is  the  range  in  which  the  curve  for  a  particle  whose  size 
makes  aQ  equal  to  10  m  The  shot  is  made  and  the  detector  outputs, 
as  photographed  by  the  oscilloscope  cameras,  are  studied  and  the  ex¬ 
perimental  distance  versus  time  values  compared  with  the  theoretical 
values.  It  can  then  be  determined  if  the  particle  was  in  the  desired 
size  range  or  if  it  was  larger.  If  it  was  in  the  desired  size  range, 
closer  comparison  with  the  theoretical  curves  will  disclose  the  actual 
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size.  A  target  placed  in  the  spray  particle  path  can  be  examined  to 
see  that  only  one  particle  was  observed. 

A  number  of  refinements  to  the  detection  equipment  will  be  suggested 
here.  First,  the  detector  box  should  be  made  larger  and  more  photo¬ 
multipliers  added.  Two  photomultipliers  are  needed  very  close  to  the 
target  area  to  measure  initial  velocity  accurately.  Another  photo¬ 
multiplier  in  the  size  measuring  range  would  give  more  accuracy  in 
plotting  the  trajectory  curve. 

Another  refinement  for  the  detection  system  would  be  to  install  a 
lens  system  in  front  of  each  photomultiplier  tube  to  collect  more  light 
and  focus  it  on  the  most  sensitive  spot  of  the  tube.  The  lens  system 
could  consist  of  one  simple  lens  with  the  proper  focal  length  and  could 
increase  the  amount  of  light  reaching  the  detector  by  almost  100  times. 

This  light  gathering  system  may  be  necessary  to  successfully  detect  one 
particle,  since  its  intensity  will  be  very  low. 

Studies  should  be  made  of  the  effects  on  spray  particle  production 

of  varying  the  impact  angle.  This  might  'be -dene-  by  rofcatfrtg~tfac  taiget - 

or  by  positioning  the  detection  equipment  at  various  angles  to  the  face 
of  the  target. 

If  k  (See  Eq.  10,  Section  II)  can  be  measured  separately  or  by 
careful  time-distance  measurements  of  a  spray  particle,  the  luminous 
efficiency  of  Eq.  17  can  be  determined.  This  would  be  of  great 
importance  to  meteor  physics. 

In  the  event  that  experimental  difficulties  prove  too  great  to 
obtain  a  single  particle,  a  cloud  of  spray  particles  could  be  studied  to 
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yield  much  useful  information  about  meteors.  A  spectrographic  study  of 
the  light  emitted  would  be  of  value. 

IV.  CONCLUSIONS 

Equations  for  deceleration,  size,  time  rate  of  change  of  mass, 
and  luminosity  can  be  written  in  terms  of  experimentally  measurable 
quantities  for  a  high-velocity  particle  moving  in  a  controlled  atmosphere. 
A  digital  computer  can  be  used  effectively  to  obtain  solutions  to  the 
equations . 

Many  important  observations  can  be  made  from  a  consideration  of 
the  solutions  to  the  equations.  It  appears  that  it  is  possible  to 
determine  the  size  of  spray  particles  by  experimentally  measuring  their 
time  of  flight  to  given  fixed  distances.  If  the  size  of  a  particle  is 
known,  its  luminosity  can  be  calculated  and  compared  with  experimental 
results  to  determine  the  validity  of  the  assumed  luminosity  equation.  If 
the  equation  is  found  to  hold,  a  value  for  t,  the  luminous  efficiency  can 
be  determined. 

Experimental  problems  may  prove  to  be  great.  It  will  be  difficult 
to  obtain  a  single  spray  particle,  and  if  it  is  obtained,  its  luminous 
intensity  may  be  so  low  that  detection  will  be  very  difficult;  however, 
detection  should  be  possible  for  many  cases  of  interest. 

Spray  particles  provide  one  of  the  best,  if  not  the  only,  means 
of  studying  meteor  physics  in  the  laboratory. 
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APPENDIX  I 


A  SOLUTION  OF  THE  EQUATIONS  OF  MOTION 
Equation  (11)  is  an  expression  for  the  deceleration  of  an  ablating 
particle. 


dv 

dt 
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— - r  v  e 


dv 

dt 


can  be  written  as 


so  that 
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Let  u  =  kv  and 
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Substituting  the  value  for  u  and  combining  equations  (a)  and  (c) 


(b) 


(c) 


du  .  3  V  '(U  '  Uo) 

ds  4  r  5  ue 


(d) 


Rearranging 


u  -  u 
e  0 


3  ^DP 

du  =  -  7  ds 

4  r  5 

o 


(e) 
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Equation  (e)  can  be  integrated  to  obtain  distance  as  a  function  of 
velocity. 


du  =  -  — 


3  CDP 


,  =  (s  -  s  ) 

4  r  o  o 

o 


(f) 


where  sq  is  the  initial  displacement  at  v  =  vq.  The  integral  in 
Eq.  (f)  may  be  separated  into  two  parts  as  follows: 


3 

4 


'(s  -so) 


(g) 


Values  are  tabulated  for  exponential  integrals  of  the  type  appearing  on 
the  left  side  of  Eq.  (g).'*’  Using  the  conventional  symbol  for  the  ex¬ 
ponential  integral  and  substituting  for  u  and  uq ,  Eq.  (g)  becpmes 
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Similar  equations  were  derived  by  Hoppe  and  Hansen  for  a  meteor 
moving  in  an  atmosphere  of  exponentially  increasing  density. 


-  50  - 


APPENDIX  II 


COMPUTER  PROGRAM  FOR  SOLVING  THE  EQUATIONS  OF  MOTION 

The  Datatron  205  digital  computer  was  programmed  to  solve  the 
following  set  of  linear  first  order  differential  equations  for  velocity 
and  displacement. 

dv  3  (~'DP  ,,2  k  (v  2  .  v2) 

TZ  ~  ~  — r  V  e  o 

dt  8  r  o 
o 


ds 

dt 


v 


with  the  initial  conditions,  s  =  t  =0. 

’  o  o 

The  program  also  provides  for  the  solution  to  the  following  equations: 
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Values  for  the  various  equation  parameters  were  put  on  cards  as 
shown . 


51 


A  new  card  can  be  punched  when  solutions  are  desired  for  any  particular 
set  of  parameters. 

The  program,  sub-routines,  and  parameter  values  are  stored  in  the 
memory  as  follows: 


Memory  Locations 


0000-0198 

0360-0479 

1010-1019 

2000-2109 

6000-6109 

The  program  is  coded  as 


differential  equation  sub-routine 

program 

parameters 

eX  sub-routine  (2.20.120) 
data  for  diff.  eq.  sub-routine 
follows : 


(2.50.010) 


MAIN  BAND 

Start  with  30  0390 


INSTRUCTION 


0360 

7360 

0  0000  34 

1000 

1 

1 

64 

5001 

2 

2 

02 

4002 

3 

3 

64 

5001 

4 

S  4 

82 

5001 

0365 

/'365 

02 

6005 

-  52  - 


0366 

7366 

0 

0000 

64 

4015 

7 

7 

82 

4015 

8 

8 

81 

6005 

9 

9 

82 

4010 

0370 

7370 

20 

0371 

1 

1 

31 

2000 

2 

2 

30 

0373 

3 

3 

82 

6005 

4 

4 

02 

6006 

0375 

7375 

64 

7019 

6 

6 

82 

4011 

7 

7 

82 

4012 

8 

8 

30 

0380 

9 

9 

1 

5037 

50 

0000 

0380 

7380 

0 

0000 

83 

4013 

1 

1 

83 

4014 

2 

2 

82 

6006 

3 

3 

02 

4001 

4 

4 

64 

7019 

0385 

7385 

81 

5001 

6 

6 

05 

7009 

7 

7 

64 

5000 

8 

8 

02 

6003 

9 

9 

20 

0171 

0390 

7390 

0 

9990 

44 

1010 

1 

1 

0 

9980 

54 

8010 

2 

2 

0 

0000 

64 

4015 

3 

3  - 

02 

0201 

4 

4 

02 

0200 

0395 

7395 

02 

0202 

6 

6 

64 

6019 

7 

7 

02 

6003 

8 

8 

30 

0440 

9 

9 

0 

5150 

00 

0000 

0400 

7400 

0 

0000 

24 

0420 

1 

1 

02 
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